The Method of Comparison Equations for Schwarzschild Black Holes 
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We employ the method of comparison equations to study the propagation of a massless minimaUy 
coupled scalar field on the Schwarzschild background. In particular, we show that this method 
allows us to obtain explicit approximate expressions for the radial modes with energy below the 
peak of the effective potential which are fairly accurate over the whole region outside the horizon. 
This case can be of particular interesting, for example, for the problem of black hole evaporation. 
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I. INTRODUCTION 

Wave equations on black hole backgrounds cannot be 
solved exactly in terms of known simple functions (for a 
review, see, e.g., Ref. Although it is usually fairly 
easy to find asymptotic expressions near the (outer) hori- 
zon or at spatial infinity, numerical analysis is what one 
eventually employs to obtain a detailed profile of the 
wave- functions. However, our understanding of the field 
dynamics on curved backgrounds would still benefit of 
analytical approximations to the exact modes, for exam- 
ple in the study of the (renormalized) energy-momentum 
tensor of quantum fields and, in particular, of the 
Hawking radiation |3] , or of quasi-normal modes Q . 

By making use of the black hole symmetry, one can 
usually separate the wave- functions and reduce the prob- 
lem to the one task of solving a second order differential 
equation for the radial part. The latter takes the form of 
the one-dimensional Schrodinger equation for the trans- 
mission through a potential barrier 



(1) 



and, for the general case of a Kerr-Newman black hole, 
is known as the Teukolsky master equation for waves of 
spin- weight s = 0, 1 and 2 [igjSj]. 

A well known method to solve equations of the form 
is the Wentzel-Kramers-Brillouin (WKB) approxima- 
tion Q , which has in fact been widely used to study 
wave equations on black hole backgrounds (for one of the 
first attempts see Ref. and for recent reviews see, e.g., 
Ref. and References therein). One of the weaknesses 
of this method is that it becomes rather cumbersome near 
the "turning points" , that is where the potential vanishes 
and one must match different branches of approximate 
solutions. The precise location of the matching points 
greatly affects the accuracy of the method and may lead 
to errors at next-to-leading order larger that those at 
leading order (for a detailed analysis of this point in the 



context of cosmology, see Ref. [lOj). The fact that the po- 
tentials in the radial equations on black hole backgrounds 
can in general have two turning points makes therefore 
it difficult to apply the WKB approximation. 

There are however improved versions of the WKB ap- 
proximation and, in the present paper, we shall con- 
centrate on the one which seems to give the most ac- 
curate leading order expressions. The method of com- 
parison equations (MCE) was independently introduced 
in Refs. [TlL IT^ and later applied, for example, to 
wave mechanics 0| and cosmology 14]. Its connec- 
tion with the Ermakov-Pinney equation was also studied 
in Ref. Its strength consists in yielding approxi- 

mate modes which are rather accurate over the entire 
domain of definition of the wave-functions. The approxi- 
mate wave-functions are however given in terms of a new 
independent variable and computing their explicit depen- 
dence on the original variable is now the most difficult 
task which requires a good deal of trial-and-crror work. 

In the next Section, we shall briefly review the equa- 
tion of motion for massless scalar fields on a spheri- 
cally symmetric black hole background. This includes 
the simplest case of linear scalar perturbations of a 
Schwarzschild black hole which will be used to show the 
application of our method in all details (including the 
use of a convenient radial coordinate and rescaling of 
the wave- function) . In Section IIIII we shall then apply 
the MCE and obtain approximate expressions for scalar 
wave-modes with energy below the peak of the poten- 
tial. This case is of particular interest, for example, for 
evaluating the grey-body factors involved in the Hawking 
effect 3]. It also presents a major complication with re- 
spect to the typical cases treated in Ref. , namely the 
presence of two zeros of the potential ("turning points"). 
We shall use units with c = G = 1 . 



II. SCALAR FIELD ON SCHWARZSCHILD 
BACKGROUND 
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Let US begin by recalling that the metric in a static and 
spherically symmetric vacuum space-time can be written 



2 



in general as 



ds^ = h{r) 



W) 



(2) 



where dO is the area element on the unit two-sphere. The 
propagation of massless, minimally coupled scalar parti- 
cles in this background is governed by the Klein-Gordon 
equation = 0. Due to the symmetry of the back- 

ground, the field $ can be decomposed into eigenmodes 
of normal frequency u and angular momentum numbers 
m as 



(3) 



where are spherical harmonics. The Klein-Gordon 
equation then separates and the dynamics is described 
by the function R which satisfies the radial equation 



_d 

dr 



■h{r) 



dR{r) 



dr 



h{r) 



R{r) = , (4) 



where is the separation constant. No exact solu- 

tion of the above equation is known even for as simple a 



background as the four-dimensional Schwarzschild space- 
time with 



h{r) — \ — f 



(5) 



where we have introduced the dimensionless radial coor- 
dinate f = t/tyv and th is the Schwarzschild radius. 

It is nowjx)nvcnicnt to introduce the "tortoise-like" 
coordinate 1241 



dx — 



f h(f) 



(6a) 



and the new radial function 



R[f{x)] = exp / h{x')dx' 



XW)\ , (6b) 



so that the radial equation takes the Schrodingcr form 



dx^ 



C^f'^x) 



t{l+\)h\f{x)\ 



1 dh\f{x)\ 

2 dx 



Xix) = 



(7) 



where we also defined the dimensionless energy u} = lu rn. 
For the four-dimensional Schwarzschild case we have 



x{f) — ln(f — 1) , 



(8) 



so that the new variable x ranges from a;(rH = 1) = — oo 
to x{'r = +oo) — +00 and x and x(x) are of Langer's 
form (7(1 . The Klein-Gordon equation finally becomes 



— +iJ^ix) 
dx"^ 



X{x) = , 



where the "frequency" 

oj^x) = {1 + e^f 



£{£+!) e^ e^ (2 + e^) 
1 + 4 (1 + e^)^ 



(9) 



,(10) 



is not necessarily a positive quantity. In fact, the above 
expression shows two qualitatively different behaviors de- 
pending on the values of uj and £. In particular, for a 
given angular momentum there exists a critical value 



'^c(^) , 



(11) 



such that if the energy lj > Qc there are no turning 
points, otherwise there are two, say uj^xi) — uj(x2) — 
with xi < X2. The latter case is the one we want to an- 
alyze in detail (see Appendix 1X1 for more details and the 
solid line in Fig. ^ for an example). 



III. MCE SOLUTIONS 

The MCE makes use of the exact solutions to a similar 
second-order differential equation. 



U[a) = 



(12) 



where Q is the "comparison function" . One can then 
represent exact solutions of Eq. in the form 



A- \ -1/2 



(13) 



provided the variables a and x are related by 



uj''{x) = 



da 

da; 



daV'^ d^ 



dx 



dx'^ \ dx 



da 



-1/2 



.(14) 



Of course, solving Eq. (|14(l is usually as difficult as the 
original problem Q, and the dependence of (t on a; can 
just be determined by using iterative schemes, in general 
cases 0,^3 or for specific problems 0,0. However, if 
we are able to find a comparison function Q sufficiently 
similar to uj in Eq. (|10|) . the second term in the right 
hand side of Eq. (|14|l will be negligible with respect to 
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FIG. 1: Comparison between the exact potential —uj'^{x) 
(solid line) and the approximate Morse potential — 0^[o-(a;)] 
obtained after solving Eqs. I19at and Ill9bll (dotted line) for 
(I) = 1/2 and £ — 2. The two turning points are xi ~ —3.13 
and X2 ~ 1.23. 



the first one, so that 



d£ 

da; 



e^(a) . 



(15) 



On selecting a pair of values xq and (Tq such that ctq = 
a{xo), the function (t(x) can then be approximated by a 
solution a = a{x) of the integral equation 



^±^dy= V±0'(p)dp, (16) 



where the signs must be chosen conveniently. This pro- 
cedure leads to the MCE approximation for x{^)^ 



da 



XMCE(a;) = ( ^ ) U{(t) , 



-1/2 



(17) 



which is valid in the whole range of x and including the 
turning points. 

We are now dealing with a problem of the form 
with two turning points, xi < X2- In order to implement 
the MCE, we need a "comparison frequency" with the 
same behavior. We shall use the Morse potential 



e'(CT) = Ae^'"-Be' 



D 



(18) 



The coefficients A and B are then fixed by imposing that 
the turning points of the comparison function are the 
same as those of the exact frequency, that is 



6(a-i = xi) = 9(a-2 = 2^2) = , 



(19a) 



and the coefficient D will then follow form the relation 



/ y/^Z^dy^ v/^e^dp. (19b) 

Jxi Jxi 

The system of non-linear equations ()19a|l and H19bp thus 
yields 



A = 



4^2 



TT^ (wi +W2 —2 y/Wl 

B = (wi +W2)A 
D = W1W2A, 



(20a) 

(20b) 

(20c) 



where the exact expression for ^ is given in Appendix IbI 
A sample plot of both the exact frequency and the com- 
parison function thus obtained is given in Fig. ^ 

The "comparison solution" is a linear combination of 
confluent hypergeometric functions |23 of the type iFi 
(for more details, see Ref. |l9j'). 



U{a) = C+ e' {^-'+^-) ^yA\+iVd + i 1 + 2 z /D, -2 z 

V2 2 VA / 

+C_e^(^°"-^")iFi ('i-zx/:D + z-^,l-2iv^,-2zVIe- 

V 2 2\/ A 

and cr(a;) is given implicitly by Eq. Hlfi|l with 

B 

^/W{p)dp^ v/0^-ln (2 A e'^ - J5 + 2 e2(cr)) ' ^ [2 D q-" - B + 20-" ^/U&Ha) 



(21) 



(22) 



r 



The above expression is rather complex. Hence, we just left of the smaller turning point (that is, for a; <C xi) 
consider the asymptotic relation between a and x to the 



(23a) 
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FIG. 2: Upper left graph: The function a = o{x) obtained by integrating Eq. 1161 numerically (solid line) and its approximate 
expression (I24II with ki = 1 and K2 = 3 (dashed line). Upper right graph: Comparison between the exact solution (solid 
line) and the approximate MCE solution 1171 with ct{x) evaluated numerically (dashed line) for C+ ~ 0.517 — 0.228 i and 
C- ~ 0.517 + 0.228 i. Panel (a) shows the percentage error for the same approximate function o given in the upper left graph. 
Panel (b) shows the percentage error for the approximate x in the upper right graph for x < xi and panel (c) for X2 < x. All 
plots are for u — 1/2, £ = 2 as in Fig. 



and to the right of the larger one (x ^ X2) 
B 



2VA 



a{x) ~ w (e^ + x) . (23b) 



As an approximate expression for a{x) for x < kiXi 
(with Ki > 1) we shall then use Eq. (|23a|) and the con- 
dition that (t(xi) = xi. Analogously, for x > K2X2 
(with ^2 > 1) we shall use the approximate solution of 
Eq. (|23bp to next-to-leading order for large x. Moreover, 
in the region between ki xi and K2 X2 we shall use an 
interpolating cubic function. The reason for introducing 
two new parameters ki and K2 ^ 1 is that the asymptotic 
forms in Eqs. H23a|) and Ij23bp may differ significantly 
from the correct a{x) around the turning points Xi and 
X2 and suitable values of ki and K2 usually improve the 
final result. To summarize, we have 



X + Xi\ 1 = I , Xq < X < KiXi 

D V VD 



Co + Cix + C2x'^ + C3 x^, Kixi < X < K2X2 (24) 



a; -f In 



X > K2X2 , 



with Ki and K2 ^ 1 and xq < ki xi is the value of x at 
which we wish to impose the initial conditions for x{^) 
and its derivative [23 • The explicit expressions for the 
coefficients Co, Ci, C2 and C3 are rather cumbersome 
and will be given in Appendix lUl 

In the upper left panel of Fig. [3 wc plot 17(2:) for the 
same values of Co and I used in Fig.Q] The solid line repre- 
sents the numerical solution of Eq. H16|) and the dashed 
line the approximate expression (|24|l with ki = 1 and 
K2 = 3. It is clear that these values of ki and K2 al- 
ready lead to a very good approximation for a around 
the turning points xi and X2 and that larger values of 
Ki or K2 would not improve significantly the final result. 
In the upper right panel of Fig. |21 we compare an MCE 
approximate wave-function with the exact solution (solid 
line) obtained numerically for the same case. The MCE 
solution shown by a dashed line is given by the expres- 
sion (|17() with U{(j) as in Eq. (|21() and a{x) determined 
by solving Eq. Hl()|) numerically. The coefficients C± are 
fixed by imposing that the MCE mode and its derivative 
equal the corresponding values of the chosen numerical 
solution &t X = —20. Since the two curves coincide al- 
most everywhere, it is clear that solving the approximate 
equation for 17(2;) is, to all practical extents, equiv- 
alent to solving the original equation 0. 
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(a) (b) (c) 

FIG. 3: Comparison between the exact solution (solid line) and the approximate MCE solution 11711 with (j{x) given in Eq. I24II 
with Ki = 1 and K2 = 3 (dashed line) for uj = 1/2 and £ = 2 (same as in Figs. Q and I^J- Panel (a) shows the result for x < ki xi, 
panel (b) for xi < x < it2 X2 and panel (c) for x > K2X2- 



In order to further clarify this point, in panel (a) of 
Fig. 121 we plot the relative difference (in percent) between 
the approximate expression 1)24(1 and the exact numerical 
solution of Eq. (|16|l . We can see that the error is less than 
10% almost everywhere, except near the zero of a{x). In 
the same figure, we also show the percentage error for 
the approximate solution shown in the right upper graph 
of Fig. 121 for X on the left of xi and on the right of X2- 
The latter error remains small, too, except around the 
zeros of x(^) due to the slight phase difference between 
the two expressions. 

Finally, in Fig. O the dashed lines represent the same 
MCE solution with a{x) given by its analytical approx- 
imation in Eq. 124(1 with ki = 1 and K2 — 3- This is 
the fully analytical expression obtained from the appli- 
cation of the MCE and is still remarkably good, with 
a more significant error around the larger turning point 
(x2 ~ 1.23). From the three graphs in Fig. and the 
upper left panel in Fig. |21 it clearly appears that such 
an error develops in the region between the two turning 
points, where the discrepancy between the exact (j{x) 
and the approximate expression ((24(1 is indeed larger. In 
this respect, let us note that, had we imposed initial con- 
ditions dA, X ^ X2, the larger error would have occurred 
around the smaller turning point {xi ~ —3.13). More- 
over, when we used a linear (or quadratic) interpolating 
function for kiXi < x < K2X2, the accuracy was (ex- 
pectedly) worse. Should one need even better accuracy, 
a fourth (or higher) order polynomial interpolation be- 
tween Ki xi and K2 X2 could be used instead. 

IV. CONCLUSIONS 

We have applied the method of comparison equations 
to a specific problem of scalar field dynamics on the 
Schwarzschild background. We have considered modes 
with energy lower than the peak of the effective poten- 
tial in the corresponding wave equation, so that there are 
two zeros (turning points). The comparison function has 
then been chosen of the Morse form with an argument 
which we have estimated analytically as well. This proce- 
dure was shown to be able to produced a fully analytical 



approximation of the wave- functions with good accuracy 
over the whole domain outside the Schwarzschild horizon. 

Let us now make a few remarks. Firstly, the applica- 
tion of the MCE is not as straightforward as the standard 
WKB approximation. In fact, the MCE usually requires 
solving some technical tasks specific to the problem at 
hand in order to obtain the argument a = cr(x) of the 
known function U — U{a). This aspect has already been 
analyzed in the simpler situation of cosmological pertur- 
bations for which one needs consider just one turning 
point In the present case, with two turning points, 
all expressions become more involved. For instance, the 
best approximation for a{x) we found convenient to dis- 
play in Eq. ((24(1 makes use of a cubic interpolation over 
the region containing the two turning points and is given 
in terms of two arbitrary parameters ki and K2 which 
can be fixed in order to minimize the error (and explicit 
expressions for the coefficients Co , . . . , C3 are given in 
Appendix [nj. Of course, one could employ higher order 
interpolating functions and further improve the result at 
the price of complexity if better accuracy is needed. 

Finally, our aim for the present work was mainly to 
test the effectiveness of the MCE to produce approxi- 
mate wave modes on black hole backgrounds. For this 
reason, no discussion of gauge fields, gravitational waves 
or spinors has been explicitly included. We in fact plan 
to apply the MCE also to such higher spin fields and fur- 
ther extend the method to study related problems, such 
as the quasi-normal modes ^ and grey-body factors Q of 
radiating black holes (also in the brane- world [22|, given 
the relevance of such objects in future collider experi- 
ments 0). 



APPENDIX A: TURNING POINTS 

If we define w = e^,we see that oj'^{x) = for the exact 
frequency in Eq. ((10(1 becomes the quartic polynomial 
equation 

w'^ + iw^ + bw"^ + cw + 1 ^ , (Ala) 



6 



with 



6 = 6- 
c = 4- 



1 



1 ^ + 1) 



2(1)2 

Its four roots can be written as 



r,2 



(Alb) 
(Ale) 



- :r6 



46- 16 -2c 7 



where i = 1, 



,4 and 



a = 12 + 62- 12 c 

/3 = 263- 366C + 27 0^ + 432 - 726 



23 a 



4 a- 



4a-^ 



23 



,(A2a) 



(A2b) 
(A2c) 

.(A2d) 



The coefficients ei = £2 = —£3 = —£4 = 9i = —62 = 6*3 = 
— ^4 = 1, and the values of Cj and i determine the signs 
of the roots. There can be at most two positive roots, in 
which case, say, ?i'3 < w;4 < < wi = e^^ < W2 = e^^, 
and only xi and X2 can be valid turning points. The 
critical value of uj in Eq. (|ll|l is then determined by the 
condition < wi = W2- 



APPENDIX B: EVALUATION OF ^ 



In order to evaluate the integral in Eq. (|19b|l . we start 
from the general relation (again with w ^ e^) 



ln{w) 



I oTTj \ [w - Wi){w - wi){w - W2) 

y w — W4 



^/{'Wl - W3){W4 - W2) 

X \ {w3 — wi){w4 — W2) E(z; W) — {w4 — wi){w4 — W2) F(z; W) 

Wi — W2 



where E, F and 11 are elliptic functions |2i, 



+ (W4 — Wi){w3 + W4 + Wi + U)2 + 2) n 
-2 {W4 - Wi){w3 + 1)(W2 + 1) n 

+2 (w4 — 'Wi)w3 W2 n 

I 

w = W2 yields 



z;W 

W4 — W2 
(W4 + l){wi - W2) 



{wi + 1)(W4 - W2) 

W4 K -W2) 

Wi {W4 — W2) 



■z-W 



z = ArcSin 



{w — Wl)(w4 ~ W2) 



[w - W4)(W1 - W2) 



(B2) 



{wi - W3){W4 - W2) 

X |(u;3 - wi){w4 - W3)E{W) 

+ {W4 - Wi) {W2 - W4) K{W) 
+ {-0)4 — Wi) {Ws + W4 + Wi + W2 + 2) n 



(Bl) 



Wi - W2 



w 



w 



{W3 - W4){wi - W2) 
{W3 - Wi){w4 - W2) 



(B3) 



-2(W3 + 1){W2 + 1){W4 - Wl) n 
+2 W3 W2 iw4 — Wl) n 



W4 — W2 
(W4 + l)(wi — W2) 
{Wi + 1)(W4 — W2) ' 
W4 {wi — W2) 



Wl {W4 — W2) 



■w 



w 

(B4) 



and Wi = e^' are the roots of u:^{x) as given in the in which we used E(7r/2: q) =E(g), F{TT/2;q) =K{q) and 
previous Appendix. The above expression evaluated at n(p;7r/2;q) = Il{p,q) |20|. 
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APPENDIX C: CUBIC INTERPOLATION for KiXi < x < K2X2. The above coefScients C^, with 

i — 0, ... ,3, can be easily expressed in terms of the pa- 
in Section Iml we use the analytic approximation for rameters A, B and D for the Morse potential |lHll as 
a{x) in Eq. (|24|l which involves the cubic interpolation 

a^Co + Cix + C2x^ + C3x^ , (CI) 



Co — KlX 



^2 a/D — 1^ K2^ X'l' (k2 a;2 — ki a;i) — 2 -JD In ^2 \fA^ k,\ x\ (3 K2 X2 — ki xi) 



2\/D{k2 X2 - Kl 

(k2 ^2 — Ki xi) K2 2^2 (2 Ki + K2 a;2) ~ 2 a/D Ki Xi (4 K2 X2 — nixi) + 12 \/15 In ( 2 Ki K2 a;i 2:2 



C2 



2\/D{k2 X2 " Kl Xl) 

6 \/D In (2 a/A ) (ki xi + K2 X2) + (2 \/D — 1 ) (ki^ xi'^ + ki K2 X1X2 — 2 K2 



^^X2^) 



c. = 



2 y/D{Ki Xi — K2 X2f' 

4 ^fD In (2 \^4) - (2 /D - l) (k2 a;2 - ki 2:1) 



2^/D{k2 X2 - Ki X\] 



(C2a) 
(C2b) 
(C2c) 
(C2d) 



Finally, one can express everything in terms of the origi- 
nal parameters Cj and I by making use of Eqs. (|20a|l - H20c|l . 
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